Superexchange via Cluster States: 
Calculations of Spin-Phonon Coupling Constants for CuGeOa 
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Calculations for spin-phonon coupling constants in CuGeOa are presented, applying fourth order 
superexchange perturbation theory to an extended Hubbard model. In our treatment, the interme- 
diate oxygen ligand states are described by band-like cluster states, due to the presence of significant 
0(p)-0(p) hopping. We also include the effect of the Ge ions on the O ligand states. 

We find a considerable q-dependence of the spin-phonon coupling. For the 7r-modes involved in 
the spin-Peierls transition, our results of the coupling constants are in fair agreement with the work 
of Werner, Gros, and Braden. Yet some discrepancies remain, especially concerning the coupling to 
the vibrations, which modulate the 0(2)-0(2)-Ge angle 5. Our studies of the pressure dependence 
of the magnetic coupling constants Ji and J2 suggest, that relatively large non-linear effects are 
present, even at small pressure values, probably due to the existence of a 'soft', van-der-Waals-type 
bond direction in CuGeOa. 
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I. INTRODUCTION 



In recent years, the magnitudes and origin of the mag- 
netic coupling constants of the quasi-one-dimensional 
quantum spin system CuGeOs have been subject of a va- 
riety of papersErO. In particular, the anti-fcrromagnetic 
first nearest neighbor coupling Ji is unusually small 
- believed to arise from superexchange hopping paths 
forming Cu-0(2)-Cu angles 77 close to 90°, which is the 
Goodenough-Kanamori-AndersonH (GKA) criterion for 
vanishing anti-ferromagnetic superexchange. 

Some authors have suggested that the hopping paths 
are strongly modified by the influence of the Ge atoms 
on the 0(2) wns, because of the strong covalent ionic Ge- 
0(2) bonds.El Other authorfil have obtained very strong 
variations of Ji with the angle 77; i.e., have predicted very 
large spin-phonon coupling constants in this system, in 
contradiction to estimates based on the data of the spin- 
Peierls transition.Q 

Very recently, Werner, Gros, and Braden (WGB) have 
presented a detailed evaluation of experimental data, 
based on a variety of theoretical estimates, to derive the 
spin coupling to the four phonon modes involved in the 
spin Pcierls transition.^ In particular, they find a strong 
coupling to variations of the angle 77, and a relatively 
large coupling to the Ge displacements. 

In the present paper we use superexchange pertur- 
bation theory to evaluate the spin coupling constants 
between the magnetic ions and their modulations with 
atomic displacements. Our niethod is an extension of ear- 
lier perturbative treatments.EI In particular, we include 
band-like oxygen ligand states, because there is signif- 
icant direct 0(2)-0(2) hopping in the chain. For our 
calculations, we start from a Hubbard- type Hamiltonian 
with five d-orbitals a, (3 per Cu site and p-orbitals at 
the O ligand sites. The full number of d orbitals at the 



Cu sites is considered in order to allow both for sym- 
metry induced orbital mixing and for the possibility of 
ferromagnetic exchange coupling involving excited Cu d 
levels. On the oxygen sites, we include the splitting and 
mixing of the oxygen p-levels caused by the Coulomb 
fields of the adjacent ions, in particular by Ge. 

The rest of the paper is organized as follows. In section 
II we discuss details of our model, in section III we first 
discuss our choices of model parameters and then present 
the results for the spin-phonon coupling constants and for 
the pressure dependence of J. Concluding remarks are 
given at the end. 



II. DETAILS OF THE MODEL 

The Cu atoms form a linear chain along the Cu02 rib- 
bons (see Fig. |l|), where the oxygen atoms in the Cu02 
plane are labeled 0(2), the apex oxygen atoms 0(1). 

We use a multi-band Hubbard-like Hamiltonian, where 
all electron-electron interaction takes place at the Cu- 
site, for details see Ref. ^. We take five d-orbitals per 
Cu atom (sites i, orbitals a) and 3 p-orbitals (sites or- 
bitals r) per oxygen and use the hole picture. We allow 
p-d hopping f^i^ between Cu ions and next nearest oxy- 
gen atoms. Also, oxygen-oxygen hopping tjj^i up to third 
nearest neighbors is included. The Hamiltonian reads 
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ef and a are the respective on-site energies of d and p 
orbitals. Uap and Iai3 are Coulomb and exchange inter- 
actions at the magnetic sites, respectively. In the atomic 
limit, the choice lap < leads to the first Hund's rule 
(maximum spin), the choice Uaa > Uap would lead to 
the equivalent of the second Hund's rule (maximum or- 
bital momentum). Electron-electron interactions at the 
ligand ions and between different sites are neglected. We 
will also use the approximation Uaa — Uaf3 — U and 
lap = Jh- We note that although this Hamiltonian does 
not yield the correct atomic multiplet spectrum, it is suf- 
ficiently general for our perturbative approach. In par- 
ticular, it incorporates the possibility for a competing 
ferromagnetic coupling of spins. 

The main impact of the Ge atoms on the electronic 
structure is to cause a significant crystal field splitting of 
the oxygen p atomic levels,u this effect is approximated 
by a crystal field Vj^"^ (including non-diagonal terms). 
Admixture of Ge 4s, 4p hole,states to the O p hole states 
appears to be rather smallllJ and thus is neglected here. 

Changes of hopping integrals (067) due to variations 
of inter-atomic distances away from the equilibrium 
separation are considered using relations like (067) ^ 
{ab^)Qe^^'^°~^\ where tq is the respective equilibrium dis- 
tance. 

The main difference to earlier perturbation theory cal- 
culations of spin coupling constants lies in the handling 
of the ligand ions. Usually, the ligands are treated in 
the atomic limit; i.e., only those hopping processes are 
considered, which occur along the spin exchange paths. 
Here we extend the method and include band-like ligand 
states. This extension is justified by the observation that 
0-0 hopping integrals of considerable magnitude exist in 
the Cu02 ribbons,EJ leading to 0(2) p-band widths of 
several eV. As we investigate the spin coupling for finite 
size Cu02 clusters, we actually deal with a finite number 
of 0(2) cluster states instead of the complete set of band 
states. 

Our calculations are carried out in a sequence of three 
steps. First, the oxygen cluster states are determined. 
As mentioned above, only 0(2) p-orbitals are considered 
in the one-particle model. Input are 0(2)-0(2) hopping 
integrals, bare 0(2) p-orbital energies and the crystal 
field caused by the adjacent ions. The cluster states are 
then obtained by diagonalizing the cluster Hamiltonian, 
using up to 16 0(2) atoms and periodic boundary condi- 
tions. In a second step, environmentally adapted Cu one- 



particle levels are determined. This calculation is carried 
out by coupling the atomic Cu levels (which are fivefold 
degenerate) of a single Cu atom to the 0(2) cluster and 
also to the apical 0(1) atoms. Here, the additional input 
are the p-d level separation and the values of pda and pd-K 
hopping integrals. The coupling is calculated in second 
order perturbation theory and leads to an effective on- 
site d-d hopping matrix, which is then diagonalized. We 
note that a rather similar d-level ordering (with respect 
to both d character and separation) can be obtained us- 
ing Coulomb-type ligand fields. In the third step, the 
adapted Cu d-levels are then used to calculate the spin 
coupling constants by usual fourth order superexchange 
perturbation theory, yet the Cu spins are transferred via 
the O cluster states instead of atomic ligand states. As 
the cluster (or band) states are delocalized, one might 
expect that the interaction between the magnetic sites 
decays rather slowly with distance. 

This procedure yields all magnetic coupling constants 
from the same order of perturbation theory, irrespective 
of the range of the coupling. As we deal with delocalized 
ligand states, instead of atomic ligand states, hopping 
processes are included up to infinite order, in principle. 
Yet, this approach remains limited to the interaction of 
pairs of magnetic sites. It should finally be noted that 
each hopping path yields a certain amplitude, and a bi- 
linear expression of such amplitudes leads to the value 
of the respective coupling constant .El The hopping paths 
may be divided into two classes, one where the interact- 
ing spins meet at a magnetic site (the conventional pic- 
ture of superexchange), and the other where they meet in 
a cluster state. The latter processes dominate, because 
of the presence of the many cluster states. Interference 
effects play a very important role, leading to some im- 
portant consequences: 

• They limit the range of coupling to essentially first 
and second nearest neighbors. 

• The 7r/2 destructive interference of the GKA rule 
is removed (actually shifted in angle). Thereby the 
crystal fields at the ligand sites, produced by the 
neighbor atoms (and napdulated by their motion) 
play an important roleJ3 

• Electron-electron interaction at the ligand sites will 
also modify the interference; however, there does 
not appear a tractable way to include it into our 
perturbative approach. 

Changes of atomic positions can be included easily in 
our calculations, this allows to determine the variations 
of the nearest and second nearest neighbor magnetic cou- 
pling constants (Ji and J2) with atomic displacements. 



III. CALCULATION OF SPIN-PHONON 
COUPLING CONSTANTS 
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A. Model Parameters 



actually turn out to be irrelevant for the values of the 
magnetic coupling constants. 



We have performed the calculations using cluster sizes 
of twelve and sixteen oxygen ions and have used three 
different sets of parameters (A, B, C). The values of the 
parameters lie in the range of values discussed in the 
literature, e.g., for the hopping integrals, we used val- 
ues similar to those given by MattheissEj As our main 
goal is to investigate the reasonable range of spin-phonon 
coupling constants, we have chosen the model param- 
eters to yield tho^jStandard values of Ji ~ 160 K and 
■hi Ji — a ^ O.S.ElBa'El The model parameters are given 
in Tab. |. Set A describes a model without coupling of 
the Cu d levels to 0(1). Sets B and C include the Cu d 
level coupling to the apex oxygens. For set C smaller O- 
O hopping parameters than in A and B have been used. 
The slight difference between 0(2)-0(2) hoppings along 
X and z is incorporated by the (3 parameters; i.e. by the 
distance dependence of ppa and ppir. We consider sets A 
and B as typical for CuGeOa, while set C with its much 
smaller ppa value may indicate the range of variations 
possible with still physically rather meaningful hopping 
parameters. 

The crystal field acting on the 0(2) p levels has been 
calculated using a point charge model for a shell of seven 
neighbor ions. These are O, Ge, Cu ions, where we as- 
sumed nominal ionic charges of —2, -f 4, -1-2, respectively 
and a dielectric constant £ = 6. For the calculation of the 
crystal field integrals, hydrogen-like 2p wave functions of 
the 0(2) sites were used. At equilibrium separation, the 
monopole contribution (1=0) to the crystal field matrix 
was removed (as this term is incorporated in the value of 
the atomic energy ep). The resulting crystal field split- 
ting of the 0(2) p levels is given in Tab. ||. We note that 
for atomic displacements, e.g., for the motion of the Ge 
ions, the monopole terms are included in the calculation 
of the spin-phonon coupling. 



Tab. Ill shows the Cu (P hole energies for the three 
sets of parameters A, B, and C for 16 oxygen ions in the 
cluster. In all cases, including the spin-phonon calcu- 
lations, the Cu (P ground state is predominantly of xz 
symmetry, with very little admixture of the other d or- 
bitals. The dxz orbital is always lowered considerably, 
caused by the strong pda hopping between the dxz hole 
orbital and the oxygen-p cluster states. All other d^ lev- 
els are much less affected. Since for set C the p-p hopping 
is considerably smaller than for sets A and B, the band 
width of the oxygen-p cluster states is reduced and, as 
a consequence, also the d level splitting. Experimental 
information on the positions of the excited Cu lev- 
els in CuGeOa is not known to us, a range between 1.5 
eV and 3 eV appears to be reasonable. Our calculations 
of the magnetic coupling constants also include the spin 
hopping paths involving excited Cu levels. However, 
these contributions turn out to be negligible, both for fer- 
romagnetic and anti-ferromagnetic exchange paths. As a 
consequence, the energy positions of the excited d^ levels 



B. Spin-Phonon Coupling 

The coupling constants in the equilibrium structure 
obtained with the three parameter sets A, B, and C are 
displayed in Tab. IV. For all sets, the results for Ji and 
J2 are very similar. Typically, the difference between 12 
and 16 site clusters is negligible. We have also deter- 
mined values for J3 and J4, here we find J3/J2 0.02 
and Ja/ ~ 0.1. Therefore, we find only J2/J1 to be 
unusually large. If we further enhance the 0(2) p band 
width by using larger p-p hopping integrals, for example, 
we find larger ratios J3/J1 and J4/J1, i.e. longer range 
coupling constants. 

For comparison with the work of WGB we have calcu- 
lated the spin-phonon coupling constants of the phonons 
involved in the spin-Peierls distortion. For this purpose 
we have determined the change of Ji with respect to 

(i) the motion of Cu along the crystallographic c di- 
rection (c = z, see Fig. |l|), 

(ii) the motion of Ge along 6, 

(iii) the motion of 0(2) along a, 

(iv) the 0(2) motion along b. 

Note that for (in) and (iv) the two 0(2) atoms bridging 
the Cu chains move opposite to each other. Also note 
that along the chain direction equivalent ions are dis- 
placed in alternating fashion (7r-mode). Our values given 
in Tab. 1^ include a factor 1/2, for comparison with WGB, 
who used a different normalization. 

A linear transformation leads to the angular coupling 
constants g,^ (to the change of the Cu-0(2)-Cu angle rf) 
and gs (to the change of the 0(2)-0(2)-Ge angle 5) and 
to the coupling to bond length variations gcu-o(2) ^-nd 



5Ge-0(2) (see Tab. VI). The coupling to the four Peierls- 
active phonons of symmetryQ can be obtained by an- 
other linear transformation, using eigenvectors of these 
modes as given by WGB (see Tab. |VIlD . 

As we have used finite atomic displacements to deter- 
mine the coupling constants of Tab. we have checked 
the accuracy of the linear transformation by directly in- 
serting the angular and bond length variations to our 
program code. These results, also given in Tab. VI, agree 
very well with the others. 

Tab. ^ indicates that three of the four coupling con- 
stants are of similar magnitude. These are 9o{2)^ 
and 5^(2) • three quantities depend on the change of 
both the Cu-0(2) bond length and on the 0(2)-Cu-0(2) 
angle rj, while is mainly determined by the Ge-0(2) 



bond length variation. As evident from Tab. VI, grj is 



much larger in magnitude than gs and is one of the dom- 
inant coupling constants in this system. This observation 
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also transforms to the result of the spin-phonon coupling 
constants of Tab. VII. Here, gi represents the couphng 
to the lowest-lying phonon mode, which mainly involves 
variations of the angle S. The coupling g2 is particularly 
strong, to a large extend this mode incorporates 77 angle 
distortions. 



In Tab. VIIl, some spin-phonon coupling constants for 

~ displacements are given (now also including values 
for J2, which are not accessible for 7r-modes). We observe 
a significant q dependence of the Ji coupling constant of 
order 15-30 %. 

We note that the WGB results differ considerably from 
our calculations. The differences may be best discussed 
looking at Tab. VI. While the values agree satisfac- 
torily, the gs values are larger by an order of magnitude 
in the WGB paper. The two other coupling constants 
9cu-o(2) and 5Ge-o(2) are comparable in magnitude, yet 
exhibit a different sign. 

We have also determined values for the spin-Peierls 
temperature Tsp using the same approach as employed 
by WGB. All our sets of coupling constants gi yield much 
smaller values of Tgp. 

We presume that the difference in the results for gs 
is the key to the discrepancies. As the (5-mode is the 
lowest-lying mode of the four Peierls active phonons, yet 
strongly involves bond-length and bond-angle rj conserv- 
ing displacements of the light oxygen ions, we suspect 
strong anharmonic contributions especially to gs- These 
contributions may show up as pseudo-harmonic in the 
WGB evaluation - thereby, via eigenvector sum rules of 
the harmonic problem, also effecting the other coupling 
constants. We note that our pressure results below also 
indicate the presence of non-linear effects, even at small 
pressures. 

An earlier calculation of g,, using exapt diagonalization 
techniques performed by Braden et alB has produced re- 
sults for the spin-phonon coupling which are considerably 
larger than our results. The calculation was performed 
using a cluster of two Cu ions, two oxygen ions, and two 
Ge ions. 

The biggest difference between our model and that of 
Ref. 1^ probably is their use of an on-site interaction Up 
at the ligand sites. This interaction affects all spin hop- 
ping amplitudes (in our perturbative picture), it there- 
fore reduces the interference effects leading to the GKA 
minimum. Additionally, there may exist a problem con- 
cerning the signs of the p-p hopping integrals. Using the 
model parameters of Ref. ^ we find reasonable agreement 
to our results only, if we change the signs of the p-p hop- 
ping integrals. (Note that our choice of sign is consistent 
- in the hole picture - with the assumption that the clus- 
ter states of predominantly a bonding type are lowest in 
energy.) 

We have also determined the value of the exchange 
constant splitting Ssp = | - |/( + Ji), i.e. the dif- 
ference between the nearest neighbor couplings in two ad- 
jacent cells in the dimerized phase. We find 63 p = 0.012. 
This is in good agreement with the results of Riera and 



Dobryja who report 5sp = 0.014. Note that 5sp ~ 0.012 
(together with J2/J1 = a = 0.35), is necessary to repro- 
duce the spin gap A = 2.1 meV.lia 



C. Pressure dependence of Ji and J2 



The distortions of the seven crystallographic pa- 
rameters, under pressure have been studied by various 
groups.E3~£j While the pressure data on a, b, and c are 
quite well established, there are somewhat larger discrep- 
ancies in the results for the four internal lattice parame- 
ters such as the value x for 0(2). These data also show 
large scattering, especially at low pressures. Addition- 
ally, we note a relatively large sample dependence of the 
p = crystallographic data. 

We have used two approaches to incorporate iJae ex- 
perimental data, especially those of Braden et al. ,0 into 
our calculation: 

(i) we use a linear interpolation 

(ii) we allow non-linear dependence of the experi- 
mental djata on pressure, using Murnaghan-type 
relations ,113 as previously proposed by Brauninger 
et al. for CuGeOa.EJ This is motivated by the ex- 
istence of the 'soft' b direction, where a van-der- 
Waals couphng dominates, see Fig. ^. 

Our results depend strongly on the kind of interpolation 
of the experimental pressure data. The p = gradients 
differ considerably (see Tab. IX), especially affected are 
lattice constants involving the 'soft' b direction. The re- 
sulting logarithmic derivatives of Ji and J2 with respect 
to pressure deviate by almost an order of magnitude, see 
Fig. |. Yet, even using the linear approach, the variation 
of Ji with pressure is rather non-linear. 

For we find -1.6 %/GPa for the linear inter- 

polation and -10.7 %/GPa for the non-linear one. In 
Ref. |l^, an experimental value of -8 %/GPa has been re- 
ported. For we obtain 1.1 %/GPa (linear) and 9.6 
%/GPa (non- linear). In a theoretical study using den- 
sity matrix renormalization group techniques, Raupach 
~ 15 %/GPa (under the assumption 

%/GPa and a - O.3).0 WGB report a value 



et al. obtain 

^ = -8 %/GPa and a - O.3).0 

of ~ — 6 %/GPa, however they extrapolate all cou- 
pling constants to the homogeneous deformation of the 
lattice from the four 7r-mode spin-phonon coupling con- 
stants, neglecting the q-dependence of the coupling. 

In summary, our pressure results are based on an in- 
complete experimental data base. Nevertheless, there are 
indications that there exists a considerable non-linearity 
in the pressure dependence, even at small pressures, 
probably due to the 'soft' b-direction of the crystal. 
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IV. CONCLUSION 

In this paper we present results of calculations for spin 
coupling constants in CuGeOa . In an extension of earlier 
work we have applied fourth order superexchange per- 
turbation theory to a Hubbard-type model with five d 
orbitals at the Cu sites and three p orbitals at the O lig- 
and sites. Because of the significant direct p-p hopping, 
delocalized 0(2) cluster states have been used. In this 
way, hopping processes up to infinite order have been in- 
cluded, equivalent, in principle, to the RKKY treatment 
of the interaction of spin pairs in metals. In our treat- 
ment we also incorporate the effects of the Ge sites, which 
mainly cause splittings of the 0(2) p levels and modulate 
them with atomic displacements. 

Our main results are: 

1 . The largest contributions to the spin-exchange cou- 
pling constants arises from those hopping paths, 
where doubly occupied cluster states are involved. 
Doubly occupied Cu states are less important. Fur- 
ther, only processes involving the Cu d ground- 
state level contribute significantly, all processes in- 
volving excited Cu d levels can be ignored. 

2. The range of spin coupling constants is limited to 
Ji and J2. Any longer distance coupling can be 
ignored for the range of parameters we used for 
our calculations. Increasing the oxygen band width 
(by increasing the p-p hopping integrals) leads to a 
longer range of the magnetic coupling. 

3. The spin-phonon coupling constants show consider- 
able q-dependence. The results for the 7r-modes are 
in satisfactory agreement with the values presented 
by WGB. There are, however, also significant dif- 
ferences, in particular, concerning the magnitude 
of the coupling to the angle S (0(2)-0(2)-Ge). We 
have investigated a certain range of 0(2)-0(2) hop- 
ping integrals (in limits which we believe are phys- 
ically meaningful), yet we could not find signifi- 
cantly better agreement with WGB results. 

4. Calculations of the pressure dependence of Ji and 
J2 yield very different results depending on the in- 
terpolation (either linear or non-linear) of the crys- 
tallographic data. The results using a non-linear 
interpolation agree much better with experimental 
data and other theoretical estimates of the pressure 
dependence of Ji and J2. 

Although the crystallographic pressure data exhibit 
significant scattering, in particular for the internal lat- 
tice parameters, and further appear to be sample de- 
pendent to some extend, we think that the non-linear 



pressure dependence may, indeed, be a real effect. The 
CuGeOa lattice structure exhibits a 'soft' direction along 
the b axis with van-der-Waals-type bonding between the 
sheets containing the Cu02 chains. Crystals with simi- 
larly anisotropic bonding, such as the spiral structures of 
Se or Te also show significant non-linearity in the pres- 
sure dataJl3 It would be highly desirable to have available 
improved crystallographic pressure data. 

Non-linear, i.e. anharmonic behavior of certain phonon 
modes may also cause some of the discrepancies between 
our results on spin phonon coupling and those of WGB. 
A candidate for large anharmonicities is the lowest lying 
mode of the four phonon modes involved in the spin- 
Peierls transition. This mode mainly consists of 0(2) 
displacements, which modulate the 0(2)-0(2)-Ge angle 
S. 

Finally, we should remark that some of the assump- 
tions underlying our studies may be questionable. This 
does not so much concern the values of the model param- 
eters - we think, we have studied a physically meaningful 
range. It concerns more the assumption that on-site in- 
teractions at the ligand sites can be neglected, since we 
find that doubly occupied cluster states are very impor- 
tant in the spin exchange paths. But note that there is 
also a dilution effect, as the extended nature of the clus- 
ter states reduces the impact of the on-site interactions. 
Unfortunately, this could not be quantified, since we do 
not see any simple way to include electronic interactions 
at ligand sites into our perturbative (cluster) treatment. 
Moreover, the magnitude of the interaction parameters 
are much less established than those for Cu^+. Exact 
diagonalization studies are limited to a few Cu and 0(2) 
atoms and also have to rely on estimates of the 0(2) 
interaction parameters. Density functional theory calcu- 
lations of Ji and J2 are hampered by the problem that 
the anti-ferromagnetic state is not the DFT ground state 
of the solid. In our view, the only way to improve the 
present level of quality in magnetic coupling calculations 
is to use ab initio quantum chemistry methods for small 
[Cu02]n clusters, properly embedded. 

In view of the above discussion we think that our 
method is a reliable tool for the calculation of magnetic 
coupling constants and related quantities such as spin- 
phonon couplings in relatively complex materials like 
CuGeOa. The main advantage is that band- like prop- 
erties of the system under consideration are taken into 
account. In addition, it allows to investigate the distance 
dependence of coupling constants in an easy manner. 
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FIG. 1. Cu02 spin chains (along z direction) with neigh- 
boring Ge atoms, x direction lies in the Cu02 plane, y per- 
pendicular to the plane. Also shown are the crystallographic 
axes a, h, and c = z. Hopping paths are denoted by the 
heavy, broken line for O-O hopping and the light, full line for 
Cu-0 hopping. In addition, the broken line between O and 
Ge indicates the ligand field of the Ge ions. 
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FIG. 2. Dependence of J\ on the Cu-0(2)-Cu angle r] for 
out of phase {q^ = tt, full line) and in phase {q^ = 0, broken 
line) distortions of neighboring cells. At the minima around 
84° (88°) Ji is slightly negative, due to interference effects 
discussed in Ref. |^. The vertical line at ~ 99° indicates the 
equilibrium angle. The dot-dashed curve shows the results of 
Ref. |. 



FIG. 3. Dependence of the 0(2)-0(2)-Ge angle 5 on pres- 
sure. The symbols show the experimental values, the full (dot- 
ted) line shows 5{p) calculated from a lineas-(|Murnaghan-like) 
fit to the experimental crystal parameters.t3 
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FIG. 4. Dependence of Ji (upper part) and J2 (lower part) 
on pressure. The squares (triangles) show the values of Ji 
( J2) calculated from the experimental structure parameters at 
the respective pressures. The full line shows Ji{p) calculated 
from a linear fit to the experimental crystal data of Ref. |l^, 
the dashed lines represent the Ji{p) of the non-linear fit. 



TABLE I. Parameter sets A, B, and C for CuGe203 in 
eV and eV/A for /3 parameters. Subscript 1 represents 0(1), 
Subscript 2 represents planar 0(2). 



u 


Jh 


ei 


pdai 


pd'Ki 


pd(T2 


pdn2 


ppa 


ppn 


/3r /3r /3r 




A 8 


-0.8 








-1.2 


0.4 


-1.0 


0.4 


1.5 2 2 


2.5 


B 8 


-0.8 


3 


-0.35 


0.4 


-1.2 


0.4 


-1.0 


0.4 


1.5 2 2 


2.5 


C 8 


-0.8 


3 


-0.08 


0.4 


-1.2 


0.4 


-0.6 


0.4 


1.5 2 2 


2.5 



TABLE IL Non-zero Crystal field matrix elements for 0(2) 
p orbitals (in eV) from the point charge model described in 
Sec.III. 

"yxx yyy yzz yxy 

-0.3 +0.2 +0.1 -0.2 



TABLE III. Energies of the symmetry adapted Cu d or- 
bitals (in eV). 





4 


4 




4 


4 


A16 


-3.5 


-0.4 


-0.1 


-0.1 


-0.1 


B16 


-3.6 


-0.8 


-0.5 


-0.4 


-0.1 


C16 


-2.4 


-0.8 


-0.5 


-0.4 


-0.1 



TABLE IV. Nearest (Ji) and next nearest (J2) neighbor 
exchange constants for CuGeOa (in K) for the sets A, B, and 
C and for 0(2) cluster sizes of 12 and 16 oxygen atoms, re- 
spectively. In addition, the results using the zero pressure 
structural parameters of Ref. [13 are given (B16p). 





yll2 


A16 


B16 


C16 


B16p 


Jl 


163 


168 


159 


158 


174 


J2 


53 


52 


49 


48 


49 



TABLE V. Spin phonon coupling constants for different 
lattice distortions (in K/A). The superscript labels the crys- 
tallographic directions a, b, c, the subscript labels the respec- 
tive ion. For comparison, the result of WGB are also shown. 



9' 


A12 


A16 


B16 


C16 


WGB 


9cu 


-271 


-278 


-274 


-319 


-890 


9Ge 


32 


34 


32 


49 


-110 


9o(2) 


397 


403 


393 


610 


400 


9o(2) 


-290 


-295 


-280 


-452 


-91 
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TABLE VI. Coupling constants to angles and distances ob- 
tained from a linear transformation of the (except gf^ and 

ddCuj which have been calculated directly from the respec- 
tive distortion pattern) in K/A. For better comparison, the 
angular coupling constants <?,, and gs are also given in K/A, 
where we related the change of Ji with angle to the distance d 
between the two outer ions defining the angle via g^ — • 



9% 


A12 


A16 


B16 


C16 


WGB 


9ri 


400 


407 


395 


566 


708 






413 


406 


582 




9i 


-15 


-16 


-0.2 


-5 


224 


9dCu 


-48 


-48 


-43 


-114 


180 


9dCu 


-48 


-47 


-47 


-116 




9dGe 


37 


38 


39 


54 


-96 



TABLE VII. Coupling constants g„,n = 1,2,3,4 {in K) of 
the four Peierls-active modes. 



9i AVI 


A16 


BIQ 


C16 


WGB 


91 -1-5 


-1.5 


-2.2 


-1.6 


-15 


92 26 


26 


26 


35.9 


58 


53 -8 


-8 


-8 


-8.2 


-30 


54 -22 


-22 


-21 


-34.1 


-12 



TABLE VIII. Spin phonon coupling constants of Ji and J2 
(in K/A) for the (qz — 0) displacements, using set B16. 





9o{2) 


9o{2) 


9Ge 


5r) 


Jl 


258 


-180 


56 


460 


J2 


-167 


119 


2 


92 



TABLE IX. Pressure gradients of the seven crystal param- 
eters at p=0 from linear and non-linear fits to the data (in 
A/GPa for a, b, and c and in 1/GPa for the internal lattice 
parameters). 





a 


b 


c 


xGe/a 


xOl/a 


x02/a 


y02/b 


lin 


-0.013 


-0.12 


-0.0045 


-0.0054 


-0.0053 


-0.0037 


-0.0015 


n-lin 


-0.021 


-0.17 


-0.0044 


-0.0093 


-0.0062 


-0.010 


-0.0024 
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